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REDUCTION THEORY AND EXPERIMENTAL 
INJECTED POLYMER SOLUTION UNDER 


INVESTIGATION OF IN SLIT 
HIGH REYNOLD NUMBERS 


Wang Xiliang Xia Changsheng 


The addition of a small amount of drag-reducing polymer 
in near-wall flow may greatly reduce the sliding friction. In 
engineering we are interested in the theory and experimentation 
of high Reynold number slit injection of drag-reducing polymer. 
Our results are presented below. 

1. Differential relations of wall shear stress of slit- 
injected drag-reducing polymer solution in flow around a stream¬ 
lined body of revolution. 


This method is based on the velocity distribution of the 
joint boundary layer in references (4) and (5). The equation 

is 


M + « «- /liny* + B + bay* + AS 

V* 


( 1 ) 


where u is the velocity component along x direction In the 

boundary layer; y+- p m y/v; p, - VtTTp ls the shear 

velocity, t is the wall shear stress, p is the density; A, 

™ v dp 

B are constants, with A = 2.5 and B = 5.5; *“■ - ,* 

dx 

v is the kinematic viscosity, p is the pressure; 

AS — FinT — r(<■) is related to the concentration 

c, v « cr is the shear velocity at the onset of brag-reduction.:. 



for zero or negative pressure 
gradient 

for positive pressure gradient. 
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The continuity equation of the flow around a body of 
revolution is 


Q(r») , d(rt') — 

6x dy (2) 


When normal stress is neglected, the equation of motion is 


u d±. + v dlL„-±.d£_ + ±QLlll 

dx dy P dx pr dy 


(3) 


where x and y are respectively the tangential and normal direc¬ 
tion of the object profile in the meridian plane; r = r Q + 
ycos<|), 4) is the angle between the x direction and the axis, 

r is the distance from a point on the object surface to the 
o 

axis. (Figure 1). 



Figure 1. Coordinates of the body of 


We assume 
with thickness 


the boundary layer to be a 

6 <<r . 
o 


Integrating equation (2), we get 


revolution. 
thin boundary layer 










[v,(A\ny* + bay* + B ■+• AZJ)) “ (.u* + A + bay*') 

dx dx 

r. * <*« , </cab)i . 


dx dx 


Substituting equations ( h ) and (5) into (3) and noticing 


that 


i<* V <*y / 


then we get 


+> ^ , * [ + y + _ J- y «-> - \-^L. v \y* (u* — ~ A ')^ 

U V * if + b *17 l y 2 ^ Jy^J t dx * * \ 2 ' 4y 

+du*\d(&B) _1 -x. J_ Kill - r; 4- ^ 

+ 'l{**-y*j?)- 4 T- -7dx + ,r dy U dx + ^dy* 


Integrating the above equation from 0 to 6 = v # 6/v we get 

iL* G + b „l *!L h _ -L i!L v \s* \(a+ ^\ b'a'S*' 

dx dx r dx !■' 2/3 

- i- AboS* - 2^1 + v% 6* («r - baS* - 2 

- vJ T s *-T r - (6) 

dx h 

where U is the flow velocity at the outer edge of the boundary 
layer; y is the viscosity; 

♦I - V 

o — « . •* — 

C - j- 0* [(<r - A)' + A' - ^ (<r - j ^)] 

H - - ^, + 2 j| V 

Let the stiffness approach infinity in equation (6) and let 
V •- V/Vm, x — x/L, r — r/L , where is the flow velocity at 

infinity. L is the length of the body. Finally we obtain the 
differential relation for the shearing stress of the body of 
revolution with slit injection. 
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(G - 3ba}» a' + £ o^B* - C) - 

+ j aS* [(/< + a - ~ b* a 1 8*‘ - -|- Abas'* - 2/i'] 

- aS* (a - 2A — baS*) (AB)' - R L V 

(7) 

in which the * and the " represent the first and second order 
derivative wrt x; 

R t - UmL/V ; 


U»)'-£ont> + r(£-£) 


2. Discussion on the differential relation of the 
shearing stress. 

1) when B=0, equation (7) is the shearing stress differ¬ 
ential relation for the flow around the body of revolution 
without the additive. 

2) From equation ( 8 ), we have (ab— T a'i/o, in plane 

uniform solution flow. Equation (7) becomes 

S*Ua- Ay + A' 4- r (a - 2A)] a - R L V or s* l(a - A )* 

+ ^-4-r (<7 - 2^)]-g-- 1 


Equation (9) has the solution 

R, - (o-„rS [<*(£+1) + < r + 2/1 > a (I+i) ~ 2A( - r ~ A ’> a i] ( 1 

when T/A is an integer. f 

• r • p 

Here a, - At?, a t _ - /< (<m e l - ) 

When r/A = 0, this is then the solution for the flow around 
the body without the additional plane. q 1 plane di 1 n f 
coefficient (2) is Cf _ J_[' L iE* _ JL [JL j* + „♦ (<r _ .♦) 4 ) 


4 


For our model, we have 


r - 2 A) (H) 

R l a 

Figure 2 shows the relational curves of C ^ to the Reynold 
number Re for poly(ethylene oxide) solution of various concen¬ 
trations T/A = 1, 2, 5 and for pure water T /A = 0. Here 
= 10 m/s, V* = 0.023 m/s. 

3) For the case of pure water, the separation condition 
for the flow around a body is G-3baH=0; in the presence of drag 
reducing polymer, the separation point condition is 

G — Iba H + rs* (<r — 2A — bA S + ) - 0 

The extra third term is usually greater than zero. When separa 
tion exists, we ^uy interpret this as moving the separation 



Figure 2. Curve of plane viscosity for uniform solution 

point further to the rear (o,6 + , G and H also change with 
additive). In practice, the curve for wall shear stress 
coefficient r m fpV l — 1/tr 1 is calculated from equation (7). 
In the case of pure water, the tail of the body of revolution 
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curve decrease since 


2 

makes o' increase rnnidly and "i /p’J 

— w 2 
V*/V and r’/r are both larr.e negative numbers. When T v ,/pU-*0, 

separation occurs. With additive, because of the extra term 

ra + (<r — 2A — bad*) s the coefficient of o' increases so that 

the fall-off phenomenon of the curve due to V* /V and r'/r is 

moderated, delaying the occurence of separation. 

4) It is very important in solving equation (7) to estab¬ 
lish the wall concentration relation along the longitudinal 
direction and to derive from which the expression for (AD)’. 
When we assume that the diffusion of the slit injected polymer 
at the boundary layer is similar to the diffusion law of the 
line source at the turbulent boundary layer, the final stage 
wall concentration 


c m GO - f/0.55 su ( 12 ) 

r 71 

where g is the unit length injection rate. When f is 

taken to be directly proportional to C w , equation (12) may be 
re-written as 


Tag* = constant 


(13) 


From (13) with the consideration of the expressions for a and 
AB, we can derive 


where 


(AB)' - [K, (<r + r 4- AB - 3 bad") - (r + AB)) 

c 

+ ir - k, cr - ib* s + )j £ - k, bas+ —■ 


K, ~ 


A B 


A B — A — baB* 


(14) 


3. Experiment on drag-reduction and wall concentration 
measurement for slit injection on a plane and on a body of 
revolution. 
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1) The measuring procedure 

The model used In this experiment consists of a plane 
and a body of revolution. The length of the plane is 3. h m. 
and its thickness is 25 rrun, with both ends trimmed sharp. There 
is a slit of width 0.8 mm at an angle of about 30° with the 
surface from the front end of the plane. During the experiment 
the part of the plane under the water line is 0.338 m. -’ix 
sampling slits of width 0.25 mm and effective length 30 - 35 
mm are placed along the longitudinal direction. Their posi¬ 
tions are tabulated in tables 1 and 2. The body of revolution 
has a length of 5.463 m. There is an annular injection slit 
of width about 0.8 mm and at an angle of about 7° with the 
surface at a distance of 0.439 m from the front end of the 
body. Seven sampling slits similar to those for the plane are 
placed along the longitudinal direction. We used the vacuum 

a 

suction sampling method with a sampling rate of about 0.3 cm. / 
s.crn. For concentration determination we used the turpidity 
method: 1 cc 50% sulphuric acid and 1 cc 1" phosphorous moly¬ 

bdenum acid are added to 10 cc of the solution and the mixture 
is shaken until the mixing is uniform. 10 minutes is allowed 
for the solution to settle and then its optical density is 
measured with a Model 72 photometer. The drag-reducing agent 
is poly(ethylene oxide) with a molecular weight of about 3 
million. The injection concentration is approximately 500 ppm 
and the injection rates for the plane and the body of revol¬ 
ution respectively are 1 kg/s and 0.74 kg/s. 

2. Experimental results 

Figures 3 and 4 represent respectively the experimental 
results of the slit injection for the plane and for the body 
of revolution. The experimental Reynold numbers are respectively 
3 x 10^ - 2.7 x 10^ and 5.5 x 10^ - 5.4 x 10^. The largest 
reduction of drag coefficient is approximately 45£ and 39%• 

For constant injection rate and injection concentration, the 
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Figure 3* Drag experiment Figure . Drag experiment 

curves of the plane for pure curve of a body c-f revolution 

water and for slit injection for slit injection of poly 

of 500 ppm poly(ethylene (ethylene oxide) solution and 

oxide). the comparison between theo¬ 

retically calculated percentag 
drag reduction and the experi¬ 
mental results. 


drag reduction decreases slightly with increasing velocity. 

This decrease is more pronounced for the plane than for the 
body of revolution because at high velocity there is less water 
immergence for the plane due to the free surface effect. 


The measured results of the wall concentration for the 

plane and for the body of revolution with slit injection are 

tabulated in Table 1 and Table 2 in which V is the flow speed 

(m/s), Q. is the injection rate (kg/s), C is the injection 
J J 

concentration (ppm), C w is the wall concentration (ppm), and 

x is the distance (m) from the front end. * indicates a higher 

concentration because the sampling bottle was not thoroughly 

cleaned after sampling for the Cj = 1000 ppm case. 

The ^results for the wall concentration shows that the 
variation of wall concentration is larger near the slit outlet. 
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The wall concentration within 1.5 n; fro::. the .oil 1 is bi,:-:er 
than the saturated concentration of 25 ppm for drar redact io: 
There seers to be a re'ion with relatively so all variation 
in the concentration for a fairly lenr distance toward the 
rear, especially for the case of the body of revolution. It 
is difficult to reduce the wall concentration to an express¬ 
ion used in reference [3], especially for the rc-'ion where 
C. = l. For the case as described in our paper, the wall ccn- 
centration must be affected by the injection ancle, the 
pressure and the surface curvature of the object. 

Table 1 Wall Concentration for the plane with injection 

— conditions ,. . ■■ —— ■ -- — ■ — ——— ; —— - ; ■■ ■ 1 - 


V-3.540 

V = 5.106 

V - 6.561 

V = 7.620 

V - 3.480 

£>, - 1.0 

£>,-1.0 

£>, = 1.0 

0, = 1.0 

£>,-1.0 

C, - 536 

C, = 536 

C, = 536 

C, = 536 

C, - 1000 


Table 2 Wall concentration for the body of revolution 

with injection. 


conditions 
of injectio 

vc. \ 


n V = 4.060 

V «= 5.SS9 

V ~ 6.962 

V — 4.230 

£>, = 0.735 

£>, = 0.735 

£>, » 0.735 

£>, - 0.735 

C, = 536 

C, = 588 

C, - 588 

C, - 1000 

104 

200 

233 

793 

116 

112 

100 

530 

114 

104 

110 

246 

• 12.6 

5.9 

7.7 

118 

1.1 

5.5 

6.7 

14.7 

6.1 

5.4 

6.2 

8.5 

6.3 

6.0 

6.2 

8.5 


t' - 4.055 
£>, •» 0.89 
C, - 588 


















4) Comparison between theoretical calculation and 
experl::.ent:al result 

We carried out theoretical calculation for the body 
of revolution of length 5.463 m at velocities = 4, € and 8 m/s. 
The method in reference [8] was used for potential flow cal¬ 
culation. Runre-Kutta method was used to solve quation (7). 

When we took C =m/0.555U and used enuation (1 J; ) for (AB) ’ as 
w ' 

derived from it, the value for the drag reduction was somewhat 

lower than the actual value. In this paper we used in our 

calculation the average value of the wall concentration as 

determined by the experiment. The initial value was determined 

from ° — wina'’’ + bae* + B + ab where the boundary layer 

thickness 5 «, q 059 S//(iogRes — 3.17) , s is the arc length, and 

( 9 ) 

the corresponding Reynold number is Res . In the calcula¬ 
tion, Afl — C.log -^ 2 - . VJhen the velocity was 4m/s, the 

curves for the shear stress coefficient r^/pU 1 — l/<^ are 
shown in Figure 5 for pure water and for poly(ethylene oxide) 
with injection concentration 500 ppm. 



2 

Figure 5- Theoretically calculated curves of x w /pU for 
water and for solution. ^ 

Integrating along the surface for wall shear stress for 
water and for solution as well as for their difference, we 
obtain respectively the value for drag coefficient and drag 
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reduction. We also plotted the theoretically calculated per¬ 
centage drag reduction In Figure from v;hlch one can see that 
the theoretical calculation basically agrees with experimental 
result. This method may be used to predict the calculated 
drag reduction for a body of revolution with injection. 
Strictly speaking, when applying Meyer's model to external 
flow, and in particular to external flow with uniform concen¬ 
tration, the assumed form AB — r(C„)ln -£*- requires proof 

fmu 


by direct measurement. 


The terms in equation (7) of this paper have clear physi¬ 
cal meaning and its calculation is simple. Therefore it fur¬ 
nished us a method to calculate drag reduction that basically 
agrees with experimental result. 


We are obliged to Comrades Zhou Dexiang and 'Wang Gueiqin 
for supplying data from concentration measurement. 
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Notation: r c , R,f )t A*nRr /ifcf^ denotes respectively the 


outer radius of the centrifugal nozzle outlet, the vortex flew ( - 

arm through which the fuel enters the vortex flow chamber, the * 

inlet cross-sectional area, and the geometrical characteris- 

p 

tic parameter; f.' , w, a, , <->=1-(r v /r g ) denote respectively j 

the discharge rate of the simple centrifugal nozzle, the dis- j 

charge coefficient, the outlet spray angle, the radius of air ; 

vortex, and the coefficient of effective outlet cross-section; 

o„ <Pi “ i ■" denote respectively the discharge rate I 

of a centrifugal nozzle with annular cross-sectional outlet, 3 

the discharge rate coefficient, the outlet spray angle, the 

radius of the central cylinder, and the coefficient of annular 'j 

outlet cross-section area; w, w , w m , w denote respectively 

a T* o 1 

the total velocity at an arbitrary point r in the outlet cross- 

sectional area, the axial velocity, the tangential velocity and 

the total velocity when the pressure energy is completely con- , 

verted into kinetic energy; p , p , p’ denote respectively the 

oil pressure in front of the nozzle, the static pressure at 

the point where the centrifugal nozzle with annular outlet 

r = r , and the static pressure at the point where the simple 
s 

centrifugal nozzle outlet r = r (all are residual pressures 

s 

relative to the ambient medium); \p is the corrected discharge 

coefficient of the nozzle inlet passage and the coefficient 

of flow column shape change as the fuel enters into the vortex 

chamber, y is the fuel weight, and g is the gravitational 

acceleration; the subscripts 1, 2, 3 denote respectively the 

results corresponding to the 3 different assumptions on p , 

s 

s denotes the value at the point where r = r or denotes the 

s 

value of the annular outlet cross-sectional area, the sub- 
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the 



I 


script v denotes the value at the air vortex 
superscript denotes the vniue corresponding to the case when 1 

the central cylinder is absent (simple centrifugal nozzle). | 


1. Centrifugal nozzels are widely used in power, aero¬ 
nautic, chemical and metallurgical industries. The operational 
characteristics of the simple centrifugal nozzle has Ions been , 

analysed by G. ’!. Abramovich [1 ] (we shall refer to it as A's 
theory). A large amount of experimental and theoretical mat¬ 
erials from analysis have been accumulated. Ordinarily for j 

furnace nozzles the experimental results do not deviate from 
A's theory significantly. This shows that the theory basically j 

reflects the operational characteristics of the centrifugal ( 

nozzle. Therefore it remains to be the foundation of indus- I 

trial calculations, and forms the point of departure for the 

[ 2 ] 

calculation of more complicated nozzles with structures. t 

I 

To increase the range of centrifugal nozzle flow rate 
regulation, various structures involving such forms as fuel 
return, nozzle outlet cross-sectional area regulation, or co- j 

axial placement of principal and auxilliary nozzles, etc. may 
be used [2,3] . In several of these cases, we encounter the 
centrifugal nozzle with annular outlet cross-section. Research 
on this kind of nozzles at present is still rather scanty. 

In this paper we shall analyse a nozzle which relies on a 
change in the cross-sectional area of its annular outlet to 
regulate the fuel rate (Figure 1). Its only difference from 
the simple centrifugal nozzle is in the installation of a 
conical-shaped pin rod at the center. The fue*l rate is regu¬ 
lated by moving the pin rod forward or backward so as to change 
the size of the annular cross-sectional area of the outlet 
nozzle, while keeping the fuel pressure constant. The typical 
result of experiments performed on the nozzle test platform 
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Figure 1. Centrifugal nozzle 
with adjustable outlet cross- 
sectional area. 

1. Orifice, 2. Vortex flow 
disc, 3. Centra] pin rod 
(can be moved forward or 
backward), fuel-separation 
ring. 


Figure 2. Regulation char¬ 
acteristics of centrifugal 
nozzle with adjustable outlet 
cross-sectional area. 

Rectangular tangential slots 
of dimension 1.87 x 2.595 
A * 1.575 


is shown in Figure 2. From Figure 2 we can see that when the 
pin rod is being moved forward and the blockage of the outlet 
is being gradually increased, there is an initial stage during 
which the nozzle flow rate and the spray angle are practically 
unaffected by the blockage. They only show significant change 
when the degree of blockage is greater than a certain value: 
a rapid decrease with increasing blockage. Similar experi¬ 
mental result is obtained in reference [3]. 

2. For conciseness, we shall base our analysis in this 
paper on the theory of an ideal fluid in a simple centrifugal 
nozzle. However, in order to obtain better accuracy in our 
calculation with the theoretical formulae, we shall first make 
some modifications to the original A's theory on simple cen¬ 
trifugal nozzles. ' 

In a simple centrifugal nozzle, owing to the centrifuge 
effect produced by the tangential entry of the fuel oil, the 
fuel does not flow through every part of the nozzle outlet 


| 












ci’oss-uect ionai area. Instead, an air vortex is formed near 
the axis of Uic nozzle, with the coefficient of effective 
c ross-s oe ti c na 1 area <p “ l-Oa/O’^l . The magnitude of <p 
is Jotor;:'.irou by A, the geometrical characteristic parameter 
of the nozzle at infinite stiffness: 


A 


x R r { 

4>fi 


( 1 ) 


The geometrical characteristic parameter A of the nozzle, 
the coefficient of effective outlet cross-section area <J>, the 
nozzle discharge coefficient and the spray angle a are related 
as follows: 


A 


b 



a 


2 sin 



( 2 ) 

(3) 

(4) 


We differ from the original A’s theory on only two counts: one 
is that we introduced in equation (1) the correction coeffici¬ 
ent to correct the effect introduced by the fact that the dis¬ 
charge coefficient in the tangential entrant passage is not 
equal to 1 and the fact that the fuel shape deforms as the 
fuel enters the vortex flow chamber, two is that we 

have taken into consideration the effect on the spray angle 
when the static pressure of the fuel is transformed into axial 
velocity at the nozzle outlet. Let sin(a/2) = (average tangen¬ 
tial velocity at the nozzle outlet)/(the total velocity when 
the pressure before the nozzle is completely transformed into 
kinetic energy) . 


Table 1 compares the experimental results with the results 
calculated from equations (1) - (4). The average value and 
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upper and lower error Units for 16 nozzles are 

•£**■ - 1.004 - 1.015 !?;!! 

while the ratio of the experimental results with that from the 
original A's theory are 

-±s*- - 0.88 «;«, -0.87 i?;* 

/***••• ««*■» 

Hence the modified or corrected formulae (1) - (^) agrees 
better with experiment than the original A’s theory. 

The nozzles listed in Table 1 all came from actual produc¬ 
tion and therefore may be regarded as representative. For 
the 2 nozzles used in our experiment, when ip is taken to be 
0.8, the calculated results from (1) - (A) also agreed well 
with experiment. When the central pin rod is ineffective 
(equivalent to a simple centrifugal nozzle), for the nozzle with 

3 tangential circular holes and A - 1.123, the discharge 
coefficient has a computational value of 0.^11, the experimen¬ 
tal value being 0.371, computed value for the spray angle is 
69.5°, the experimental value being 61°; For the nozzle with 

4 rectangular tangential slot and A = 1.575, the computed value 
for discharge coefficient is 0.337, the experimental value 
being 0.335, the computed value for the spray angle is 78.7°, 
the experimental value being 72°. 

In Table 1, we take = t.8 to compute A. Strictly 
speaking, ip should be dependent upon the structure of the 
nozzle and its state of mechanical processing. For simplicity, 
we shall take ip = 0.8 in all the subsequent analysis. 

For computation of equations (2) and (3) it is possible 
to use the curve in reference [1], Equation (2) - (4) may also 







-MU 


Table 1. 



Key: (1) Nozzle No.; (2) Average; (3) Geometric Character¬ 
istics; (4) Discharge rate; (5) Ratio of experimental 
to theoretical; (6) Spray angle. 


be approximated with the following relations: 


0.428 A-* M 

(5) 

0.630 

(6) 

65.2° A 0 M7 

(7) 


We obtained these approximate relations by using the least 
square method on a numerical calculator. In the region 
A=0.7 - 4.0, the results do not differ from those obtained by 
using equations (2) - (4) by more than 2%. 


3. Let us now analyse the operational characteristics 

of the centrifugal nozzle with an annular outlet. Analogous 

to the treatment in reference [1], we consider the case of 

an ideal fluid. Let the radius of the nozzle outlet be r , 

c 

and the radius of the concentric cylinder be r g so that the 

actual nozzle outlet cross-sectional area is the annular region 

bounded by the concentric circles with radii r and r (Figure 

c s 

3). There are two possible cases: (1) r is less than the 

s 

radius of the air vortex as calculated from the A value; (2) 
r g is larger than this vortex radius. We discuss the cases 




17 

































separately below: 

1) Case of r y (Figure 3a). Here the cylinder only 

occupies a part or the whole of the vortex cross-sectional 
area, and will not affect the flow of the oil. Therefore the 
discharge coefficient and the spray angle of the nozzle are 
exactly the same as in the absence of the cylinder. VJe can use 
the formulae (2) - (4) for the simple centrifugal nozzle in 
our computation. This prediction has been verified by our 
experiment: from Figure 2, when (r,/r,y<0.38 , the dis¬ 

charge coefficient and the spray angle of the nozzle are 
independent of (r /r ) 2 when the computed value of (r /r ) 2 
is 0.384, which demonstrates that for r g < r , the magnitude 
of r g does not affect the discharge coefficient and the spray 
angle of the nozzle. The experimental result of reference [3] 
agrees with this conclusion. 

2) Case of r g >r v (Figure 3b). Here the insertion of the 
central cylinder not only will make the air vortex disappear, 




Figure 3. Outlet Cross-section of a centrifugal nozzle 
with an annular outlet 
(a) r s <r v (b) r g > r y 


but will also directly affect the area of the oil flow. Let 
G and respectively represent the flew rate and the discharge 








coefficient- of the c^ntr1 t’ugai no-^lo with the annular outlet 

p be the residual pressure of the fuel ell in the nozzle 
o 

relative to the ambient .'tedium, y be the weight: of the fuel 
oil, arid £ be the gravitational acceleration, then 


C. 


h‘* r W 2 


nrp» 


( 8 ) 


It is easy to prove with the same method as in reference [1] 
that: the axial velocity in the cross-sectional area of 

the outlet is still uniformly distributed, independent of r. 

G, 

Thus, w t — -;—, where 


y uioi-i j- uu t' 


r* r.<p, 

We notice that the total velocity at the inside boundary r 

s 

of the nozzle's annualr outlet (assuming P s to be the resi¬ 
dual static pressure with respect to the ambience at that 
point) is 


in which the axial velocity 
tial velocity -- OlA - G> 



and tangen- 


Wj, — 


rfi4>r, V*r\J 1 — 


Hence 


I - Jl 

G, — r* T \ o —t ==^= j= 

V 


(9) 


Comparing (8) and (9), we know that the discharge coefficient 
of the centrifugal nozzle with annular outlet when r >r is 

S V 


»-'F-f.Nh + : t=7. 


For a given nozzle set of and operating conditions, p , A, <)> 
are all known quantities. Therefore the key to computing u g 
is the determination of the residual pressure p at the surface 

o 

r . Generally speaking, p may be a complicated function of 
s s 

<J> (or r ). Since at present there is a lack in direct experi- 
s s 

mental data about p , we adopt the method of making all kinds 

s 

of assumptions on p , deriving the corresponding formula for 

s 


T 
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U and then commaring the different commuted results with the 
experiment to determine which assumption about p is more in 
agreement with reality. 


We have made the following three assumptions about p g : 

(1) We assume that the static pressure at r=r c ^ stays the same 
before and after intreducing the central cylinder. Thus the 
static pressure p o ^ at tho surface of the cylinder is equal to 
p’, the static pressure at r=r in a simple centrifugal nozzle. 
Analogous to reference [1] 

_ i - „?.) - i(£*Y(i _ X) _ <„(_> -_J_) 

2g 2g 'Y'l’fi' ' r ti rj/ • \ 1 — V 1 — <p*' 

whence 


fa 


Ji-to-li- - _ 1 _V 

Y pt V v 1 — V 1 — 9V 

Substituting into equation (10), we get 

- t = t ) 

JT73Z 

V V? 1-9. 


( 11 ) 


(2) We assume that the axial velocity at the nozzle outlet 

stays constant before and after the introduction of the cylinder, 

i.e. w /w’ = 1. Since . . , . 

a a JZ*- - {<L\ (V\ HdSL 


w, \Gv \<p,/ n<p. 


then we have 


/‘ai *“ /‘<P«/<P 


( 12 ) 


Substituting equation (12) into (10), we see that this is 
equivalent to assuming 
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0 on the surface of the 


( 3 ) We assume that the pressure p 0 = 

central cylinder at all times, i.e. even r >r , the static 

s v 

pressure on the cylindrical surface is always equal to the 
ambient pressure. From equation (10) we knew 





(13) 


All of the equations (11) - (13) derived from the three differ¬ 
ent assumptions give results in qualitative agreement with the 

experiment, namely that when r >r , as r is increased (d> is 

s V s s 

decreased), both the discharge coefficient and the spray angle 

of the centrifugal nozzle with annular outlet decrease in value 

(see below). To determine which assumption about p is in 

s 

quantitative agreement with the experiment, we compare the com¬ 
puted results of equation (11) - (13) with experiment. Figure 
4 and Figure 5 present the ratio of the computed to the experi¬ 
mental results for the tangential circular hole type nozzle and 
the tangential rectangular slot type nozzle respectively. 

From Figure 4 we can see that for this tangential hole type 

nozzle, when <f> is large (small blockage), the experimental 
s 

points follow more closely to the curve calculated from equa¬ 
tion ( 12 ) and when <f> decreases (blockage increases), the experi¬ 
mental points gradually shift over to the curve of equation ( 11 ). 
From Figure 5 we can see that for the tangential slot type 
nozzle the experimental points follow more closely the calculated 
curve of equation ( 11 ). 

The author of reference [ 3 ] changed <J> by installing con- 

s 

centric cylinders of different diameters in the center of the 
tangential hole type nozzle when he investigated the problem 
of the Interaction between the principal and auxiliary nozzles 
in a gas-fired turbine. His result showed that the experimental 


points follow the curve of equni icr. ( ] 2 ){Figure 6) . 



0.2 0.4 0.6 08 1.0 



0.2 0.4 0.6 0.8 18 

<P. 


1. Computed from (11) 2 from (12) 3 from (13) Ratio of 

computed to experimental results 


Figure ^. Discharge coeffi¬ 
cient of centrifugal nozzle 
with adjustable outlet area 
of the tangential hole type. 

three #3.2 tangential circu¬ 
lar holes A = 1.123 


Figure 5. Discharge coeffi¬ 
cient of centrifugal nozzle 
with adjustable outlet area 
of tangential slot type. 

four rectangular tangential 
slots 1.8? x 2.595 A --- 1.575 


Similar to the above, for r > r the spray angle is calculated 

S V 


from the following equation: 


*(?H 


where w^ is the average tangential velocity at the outlet of 

the centrifugal nozzle with annular outlet, which we take to 

be the tangential velocity at (r + r )/2, i.e. 

s c 


. - £JL _ _ Idfh _ lit 

1 + J 1 ~<p< ' r 


The total velocity w when all the pressure energy is trans¬ 
formed into kinetic energy is 


■ i^ki ijj i 





Substitute the different i: c , derived from the different 
assumptions (equation 11 - 13) into equation 14, we get the 
corresponding formulae for a : 


u) r s = p, s *->:■*» 


o.l 


— 2«n~* ( ■ ) 

n + Vi - <p/ 


(15) 


(2) Constant axial velocity before and after introducing 
the cylinder , . 

„ • -i / 2/ ;<,» \ 

(16> 


(3) P s "0 


o.l 


2 sin"* (- 

V l + V1 - <p/ 


(17) 



Figure 6. Ratio of discharge coefficient of centrifugal nozzle 
of tangential orifice type with annular outlet to the computed 
result from equation 12. 

four circular inlet holes 





Hence 




2 sin 


-i 7 7Ali ' _ \ 

'1 + \/l — <pj 


(14) 


Substitute the different \i^ derived from the different 

o 

assumptions (equation 11 - 13) into equation 14, we get the 

corresponding formulae for a : 

s 


(1) P s = p ’s P.-P.tt> - 2 sin-* ( ■ - - ) 

N 1 + VI — <P. 


(15) 


(2) Constant axial velocity before and after introducing 
the cylinder 

2 A f‘,i \ 


(3) P s =0 


„ — 2 sin"' (- f 

'1+vl - <P.' 

a„ - 2*in"‘ (- 

V J + VI - W 


( 16 ) 


(17) 



Figure 6. Ratio of discharge coefficient of centrifugal nozzle 
of tangential orifice type with annular outlet to the computed 
result from equation 12. 

four circular inlet holes 
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For Uk 


ntial hole typo and the tangential slot type 


no::::leu that we used in our experiments, the ratio of’ the 
experimental result to that computed from equation (15) - (17) 


are shown 


ely in Figures 7 ant 


8. It can bi 


;een that 


the different formulae all give results in qualitative agreement 

with experiment, na::.ely that when decreases, the spray angle 

a also decreases. Fxeeot for small values of 4 when the 
s ■ s 

error is relatively large, the experimental points for the tan¬ 
gential hole type nozzle are closer to the computed result from 
equation (1.6) while those for the tangential slot type nozzle 
are closer to the computed result from equation (15) (Figure 8). 
The principal reason why the computed values deviate from 
experiment for small 8 g is that the parts of the nozzle are not 
truly concentric. The nozzles used in our experiment are used 
in actual production. There is always some degree of eccentri¬ 
city between the central pin rod and the orifice. When the 
blockage is large (or when <6 is small), the spray cone deviate 
significantly from axial symmetry so that in some directions 
the spray angle is much larger than that in some other directions. 
This also affects the discharge characteristics of the nozzle, 
and may well be what caused the discrepancy in the results for 

small <f> in Figure A and 6. 
s 


Thus, for the tangential hold type nozzle, the discharge 

coefficient (Figure A and 6) and the spray angle (Figure 7) 

both support the assumption that the introduction of the central 

cylinder does not affect the axial velocity at the nozzle outlet; 

however, for the tangential slot type nozzle, the discharge 

rate (Figure 5) and the spray angle (Figure 8) both support the 

assumption that the static pressure at r is not affected by 

s 

the introduction of the central cylinder. Further investigation 
is necessary to find out why there is a difference in the regula¬ 
tion characteristics of the two types of nozzles with different 
structures. 


2A 





1 is calculated Kith equation (15); 2 with equation (16); 
3 with equation (17) 


Figure 7. The ratio of 
calculated to experi¬ 
mental values of the spray 
angle for a centrifugal 
nozzle of the tangential 
hole typo with adjustable 
outlet cross-sectional 
area. 

Three 3.2. tangential holes 


Figure 8. The ratio of 
calculated to experimental 
values of the spray angle 
for a centrifugal nozzle of 
the tangential slot type with 
adjustable outlet cross- 
sectional area. 


Four 1.87 x 2.595 rectangular 
tangential slots. 



4. In summary, the operational characteristics of the 
centrlfural nozzle with annular outlet are ar follows: 

1) When the radius of the central cylinder is not larger 
than the value of r computed from the value of A (equation 1), 
the introduction of the central cylinder does not affect either 
the discharge rate or the spray angle of the nozzle. Computa¬ 
tions may be carried out with the equations for the simple 
centrifugal nozzle, equations (2 - 4) (or c rial:, ns (5 & 7)). 
For the nozzle with co-axially arranged principal and auxiliary 
nozzles, the outer diameter of the auxllliary orifice should no 
be larger than the air vortex of the principal nozzle (calcula¬ 
ted from equation 2). 

2) When r >r , as <J> decreases (blockage increases), the 
discharge coefficient and the spray angle of the centrifugal 
nozzle with annualr outlet both decrease. The decrease of the 
discharge rate is necessary for regulating the fuel rate, but 
the decrease in the spray angle is usually undesirable. It is 
just because of this that not too many nozzles with only adjust 
able nozzle outlet cross-sectional area are used in practice. 
For such nozzles, the range of the regulation ratio should not 
be too large, i.e. <J> should not be made too small, so that too 
small a spray angle may cause an overly concentrated fuel dis¬ 
tribution and poor atomization. Within the practical range of 
regulation ratio, we suggest that equation (12) and (36) be 
used to compute the discharge rate and the spray angle for a 
tangential hole type centrifugal nozzle with adjustable outlet 
cross-sectional area and that equation (11) and (15) be used 
for the tangential slot type nozzle. 

If the nozzle inlet cross-sectional area is also adjusted 
at the same time when the outlet cross-sectional area is 
adjusted, then the problem of the changing spray angle when 
only the outlet cross-sectional area is adjusted may be over- 
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come. If d-.-s li.n-.-J j.-rlv, the spray an.:!-- may wo 1 1 remain 
unchanged win-n *-1/' discharge rate of t he nozzle Is In* inn 
adjusted. Some exr erl mental result in this direction was given 
in reference [•'*], inJicatinr that this is really an improvement 
with practical value, h w- v- r, i--.’'-r-n.-- [hj failed to give the 
correct ec::;u\ at! onnl so heme for- this case of siniultaneously 
adjusting both the inlet and outlet cross-sectional areas. A 
slight modification in the computational scheme in our paper 
should enable us to apply it to this more complicated situa¬ 
tion; we need only to consider that in the regulating process, 
as varies, the inlet cross-sectional area also varies, hence 
the geometrical characteristic parameter A calculated from 
equation (1) is now a variable. The effect of this computational 
scheme will be left for future investigation. 

Addendum: The cool state experimental result of the power 

plant fuel nozzle cited in Table 1 was obtained jointly by the 
Beijing Central Research Center, Department of Hydroelectricity, 
the Beijing Thermal Power Generating Station and the Xian Ther¬ 
mal Engineering Research Center. We also cited part of the 
cool state experimental result on nozzles with adjustable out¬ 
let cross-sectional area done at the Capital Steel Forging 
Plant. Comrades Jin Yongli, Zhang Buzhou, Zui Xiansheng, Zhang 
Kerning, Huang Jiaqi, Yu Zhemin and Wang Shengli also partici¬ 
pated in this experiment. 
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